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ABSTRACT 

We present deep H I 21-cm and optical observations of the face-on spiral galaxy M 83 obtained as 
part of a project to search for high-velocity clouds (HVCs) in nearby galaxies. Anomalous-velocity 
neutral gas is detected toward M 83, with 5.6xl0 7 Mq of H I contained in a disk rotating 40-50 
km s _1 more slowly in projection than the bulk of the gas. We interpret this as a vertically extended 
thick disk of neutral material, containing 5.5% of the total H I within the central 8 kpc. Using an 
automated source detection algorithm to search for small-scale H I emission features, we find eight 
distinct, anomalous- velocity H I clouds with masses ranging from 7xl0 5 to 1.5xl0 7 Mq and velocities 
differing by up to 200 km s _1 compared to the H I disk. Large on-disk structures are coincident with 
the optical spiral arms, while unresolved off-disk clouds contain no diffuse optical emission down to a 
limit of 27 r' mag per square arcsec. The diversity of the thick H I disk and larger clouds suggests 
the influence of multiple formation mechanisms, with a galactic fountain responsible for the slowly- 
rotating disk and on-disk discrete clouds, and tidal effects responsible for off-disk cloud production. 
The mass and kinetic energy of the H I clouds are consistent with the mass exchange rate predicted 
by the galactic fountain model. If the HVC population in M 83 is similar to that in our own Galaxy, 
then the Galactic HVCs must be distributed within a radius of less than 25 kpc. 
Subject headings: galaxies: individual (M 83) — galaxies: ISM — galaxies: kinematics and dynamics 
— radio lines: galaxies 



1. INTRODUCTION 

Our understanding of the gaseous content of our 
Galaxy has advanced greatly in recent years as instru- 
ments sensitive to a variety of wavebands have become 
available. One component of the Galaxy's gaseous neigh- 
borhood remains poorly understood: the ensem ble of 
hi gh- velocity clouds (HVCs). Initially discovered in ll963l 
bv lMuller. Port, fe Raimondl . HVCs are clouds of H I gas 
that deviate in velocity from the differentially rotating 
disk of gas by 70-300 km s _1 . These structures have 
been studied extensively in H I emission, Ha emission, 
and absorption against background point sources, with 
several pub lished catalogs detailing their observed char- 
acteristics dWakker fc van Woerdenl ll99lL fl997t [Wakker! 
[2001: IPutman et alJl2003T ~ 

A cursory inspection of the HVC catalogs reveals a 
tremendous variation in size, morphology, and velocity; 
this variation has led to a number of formation sce- 
narios that predict different physical properties. One 
favored model is the Galactic fountain, a cyclic phe- 
nomenon whereby hot (T ~ 10 6 K) gas is ejected into 
the lower halo (z < 10 kpc) by supe rnova explosions 
((Shapiro fc Field! 119761 : iBregmanl [l980h . Here the gas 
cools and loses its buoyancy, free-falling back to the disk 
in sheets at mainly negative velocities approaching 100 
km s _1 . A second scenario suggests accretion of close 
companion galaxies, either through tidal or ram-pressure 
stripping as a companion moves through the hot, ex- 
tended Galactic halo. This is undoubtedly the origin 
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of the Magellanic Stream (|Mathewson et aLlll974h . and 
there is rece nt evidence for gaseous remnants of satel- 
lite galaxies (|LockmajJl2003t IPutman et al.ll2004f l. Other 
scenarios propose alternate extragalactic origins, such 
as ac c reted interg alactic medium (IGM) material (jOortl 
119661 I1970L fl 981 l) or discrete remnants of galaxy for- 
matio n (jOortJ 119661 : IBlitz et all 119991 : IBraun fe BurtonI 
1999). This last model attributes the most compact 
HVCs (CHVCs) to an ensemble of primordial, dark- 
matter-dominated "minihalos" streaming into the Local 
Group along filaments of moderate overdensity. Such 
clouds would have a typical distance of ~ 750 kpc, 
and total mass of ~ 10 8 -10 9 Mq per cloud (as s uming 
M H i = O.lMtot)- Modeling bv ISternberg et al.l (120021) 
and si mula tions of CDM halo evolution bv iKlypin et all 
(1999) and lMoore et alj (|1999f) support a more local dis- 
tribution of ~ 150 kpc. The clouds need not be dom- 
inated by dark matter, as uneven CDM halo cooling 
can result in fragmentation and produce a population of 
moderately-ioni zed, pressure-confined c louds embedded 
in a hot corona ((Mailer fc Bullock! 12001 . 

One key discriminator between these production mech- 
anisms is distance, yet this parameter is difficult to de- 
termine since HVCs contain no visible stars. The most 
direct method of searching for HVC absorption against 
background halo stars has met with limited success until 
very recently, yielding upper limits on the distances to 
five HVCs and eight intermediate- velocity clouds (IVCs; 
IWakkerl [200ll) . The upper bounds to the vertical z 
heights of these HVCs range from 0.2 kpc to 7 kpc, and 
the corresponding mass limits range from about 1 Mq for 
a very small HVC to 2xl0 6 Mq for Complex A. Recent 
work has placed stronger constraints on an additional 
four HVCs, placing them at distances of ~ 5-15 kpc 
(iThomet all I2006U Wakker et al.l I2007L 12008( 1 . In par- 
ticular. IWakker et all (|2008f l infer a distance of 3.7-11.2 
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TABLE 1 

Predictions of HVC Production Scenarios 



scenario 


distance 


kdevl 


size 


■Mcloud 


z 




(kpc) 


(km s" 1 ) 


(kpc) 




(Zq) 


IGM infall 


< 3 


< 150 


< 1 


<10 4 


0.1-0.3 


Galactic fountain 


< 10 


< 150 


< 1 


10 4 -10 5 


> 1 


companion stripping 


5-100 


< 300 


~ 1 


10 5 -10 6 


0.1-1 


circumgalactic warm clouds 


~ 150 


< 300 


~ 1 


10 6 -10 7 


0.1-0.3 


circumgalactic DM halos 


~ 150 


< 300 


~ 1 


10 6 -10 7 


0.1-0.3 


Local Group DM halos 


~ 750 


< 300 


1-10 


10 8 


0.1-0.3 



NOTE. — Sizes and masses shown are typical values based on the median 
observed parameters of the HVC ensemble. Masses for the first three scenarios 
assume H I is the primary constituent of the cloud, while the others assume an 
neutral gas fraction of 0.1 (warm clouds) or H I to dark matter ratio of 0.1 (dark 
matter halos). Velocities for the Galactic fountain and IGM infall models assume 
an adiabatic corona of fully ionized gas with T ~ 10 6 K. 



kpc for Complex C, suggesting an H I mass of 3-14 xlO 6 
M Q . 

Indirect methods have also been used to deter- 
mine HVC distances, such as constraining the de- 
gree of ionization on the cloud surface. It is diffi- 
cult to model the strength of the UV radiation field 
leaking out of the Milky Way, although there are 
indications that the field is sufficient to ionize at 
least t he surfaces of an HVC (iBregman fe Harrington! 
119861: iBland-Hawthorn fc Malonevl Il997ft ~ Several 
groups have i dentified ionized "skins" in a num- 
ber of HVCs dWeiner fc Williams! fl99fit iTufte et all 
199& iBland-Hawthorn et al.lll998h . and observations by 



have uncovered anomalous H I gas in every target im- 
aged with sufficient depth. Deep radio synthesis observa- 
tions of NGC 891 (ISwaters et alj|1997t lOosterloo et all 



Putman et all (|2003h of Ha emission in 25 HVCs limit 



their distances to 5 < z < 40 kpc. 

While recent progress has been made, large gaps re- 
main in our understanding of Galactic HVCs, the nature 
of the gaseous Galactic halo, and the relationship be- 
tween such a halo and the thin, star-forming disk. The 
various HVC production scenarios predict different con- 
sequences for galaxy formation and evolution. In Table 
[TJ we outline the predictions of the scenarios described 
above in terms of observables (size, velocity) and phys- 
ical properties (mass, metallicity) . Our vantage point 
remains an impediment to discriminating between these 
scenarios, as we are embedded within the very medium 
we attempt to disentangle. The superposition of mate- 
rial in position and velocity complicates the analysis, and 
chance locations of statistically anomalous sources (such 
as the Magellanic Stream and large HVC complexes) bias 
the conclusions. 

These problems are largely overcome by observing a 
large sample of external galaxies, constraining the fun- 
damental properties of their diffuse halos, and extrap- 
olating to the Milky Way. Until recently, extragalactic 
H I studies have concentrated on the structure of the 
kincmatically cold disk, and the discovery of anomalous- 
velocity material has been serendipitous. Early evi- 
dence of anomalous H I appeared in th e face-on galaxies 
M 10 1 (jvan der Hulst fc Sancis] 11988) : iKamphuis et all 
[1991 NGC 628 (IKamphuis fe Briggsl fl992l) . and NGC 
6946 (|Kamphuis fe Sancisil 119931 ). with detections of a 
few clouds of mass A/hi > 10 7 Mq. Integrated H I pro- 
files from a sample of 14 face-on galaxie s show evidence 
for hi gh- velocity gas in 10 of the systems (|Schulman et al.l 
il994f ). More recent observations of inclined galaxies 



[20071). NGC 2403 ilFraternali et al.M2002bD. NGC 4559 



(|Barbieri et aT1l2005l), NGC 253 (iBoomsma et al.ll2005h 



and NGC 6946 (jBoomsma et al.ll2008f ) have revealed ver- 
tically extended, slowly rotating "thick" disks of H I, 
suggestive of Galactic fountain activity. Modeling of a 
subset of these data confirm this conclusion, although 
some amount of IGM infall is required to prevent escape 
of the fountain gas and to produce the gas kinematics 
that are observed (jFraternali fc Binnevl 12008). A GBT 
study of M 31 has uncovered for the first time a popula- 
tion of discrete HVC -like H I clouds in an external galaxy 
(jThilker et al.ll2004D . These objects are all within 50 kpc 
of the galaxy. Surveys of galaxy groups have so far failed 
to detect HVC analogs without o ptical counterparts (e.g., 
IZwaan 2001; Pisano et al.ll2007l and references therein). 
Extrapolation by iPisano et al.l f"2007) of results from six 
Local Group analogs suggest the Galactic HVCs all lie 
within 90 kpc of the Galaxy. 

In this paper, we describe the first results of a search 
for HVCs in nearby external spiral galaxies. We concen- 
trate on face-on galaxies, since the line-of-sight velocity 
of the rotating H I disk is minimized and the parameter 
space available for HVC searching is maximized. These 
targets provide a necessary complement to the numer- 
ous inclined spirals which have been recently observed, 
as detailed above. By combining the H I data with deep 
optical surface photometry, we can constrain the amount 
of starlight in any anomalous H I clouds we find. 

The next section summarizes the observations and re- 
duction of both the H I and broad-band optical data. 
In Section 3, we discuss features of the bulk H I and 
optical disks. In Section 4, the cold H I disk is mod- 
eled and subtracted, and extended anomalous H I emis- 
sion is analyzed. Section 5 discusses the detection of 
discrete anomalous H I emission, introducing a new de- 
tection algorithm (further detailed in the Appendix [A| 
and presenting simulations that constrain the effective- 
ness of this technique. The discussion in Section 6 places 
our results in context of previous work and predictions. 
Our conclusions are summarized in Section 7. 



2. OBSERVATIONS & DATA REDUCTION 
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TABLE 2 

M 83 Properties 



NGC number 


5236 


Type c 


SBc(s)II 


RA (J2000) 


13h37m00.8s 


Inclination 11 


24° 


Dec (J2000) 


-29°51'59" 


Pos. angle d 


225° 


I 


314.6° 


D 25 h 


12.'9 


b 


32.0° 




17 kpc 
7.'9 


Distance a 


4.5 Mpc 


^Holm 


^hel 


503-516 km s" 1 




10 kpc 


B h 


8.2 


Total mass c 


1-5 xlO 11 A/ Q 


M B 


-20.1 


HI mass c 


6.1 xlO 9 M 


A v 


0.22 


1' 


1.3 kpc 



a IKarachentsev et alj H2002I') 
b Ide Vaucouleurs et at] l l 199111 
c ISandage fc Tammannl Tl987) 
d ITalbot et al.l 

c Huchtmeicr & Bohncnstcngcl (1981) 



M 83 is a nearby (4.5 Mpc; IKarachentsev et al.ll2002f ), 
face-on grand-design spiral galaxy located in the Centau- 
rus A group. Wh ile the optical, star- forming disk is small 
and well defined (|Talbot et al.l(l979fl . the H I disk is very 
extended, reaching to 6.5 times the optical Holmberg ra- 
dius with 80% of the total H I mass found outside this ra- 
dius (|Huchtmeier fc Bohnenstengel|[l98ll ). The H I disk 
is resolved into distinct rings and arms in this region, ex - 
hibiting a high degree of warping (|Tilanus fc Allen| [l993) . 
The basic properties of M 83 are presented in Table [2j 

2.1. HI 21-cm Data 

The radio data for M 83 were obtained in three sepa- 
rate observations with the VLA 4 on 24 and 28 February 
and 1 March 1999. To match the beam size with the 
scale of the expected emitting regions (1 kpc = 46" at 
4.5 Mpc), we chose to observe in the DnC configuration, 
in which the northern arm of the array is in the more 
extended C configuration. For sources at southern decli- 
nations such as M 83 (meridian altitude = 26? 1 at the 
VLA's latitude of +34?1), this configuration produces a 
more symmetric ww-plane coverage pattern and a synthe- 
sized beam of about 35" FWHM (0.76 kpc at a distance 
of 4.5 Mpc). Each observing run consisted of a 5 min 
observation of the primary flux calibrator 3C 286, with 
a scan integration time of 20 s. This was followed with 
a series of ~ 50 min exposures of M 83 with integration 
times of 60 s, intermixed with 5 min exposures of the 
secondary calibrator 3C 283. The total time spent on 
source was 12.2 h. 

The VLA spectral line correlator was configured to ob- 
serve in the 1420 MHz frequency band, or "L-band" . 
The correlator configuration and resulting data param- 
eters are summarized in Table [3l To allow detection of 
H I emission at highly anomalous velocities, we used the 
largest bandwidth offered by the correlator, 3.125 MHz. 
This corresponds to a velocity width of 660 km s _1 , of 
which only the inner 540 km s _1 were retained because 
of edge noise. With the bandwidth centered on the sys- 
temic velocity of M 83, this enabled detection of fea- 

4 The VLA is operated by the National Radio Astronomy Ob- 
servatory, a facility of the National Science Foundation operated 
under cooperative agreement by Associated Universities, Inc. 



tures deviating by ±270 km s _1 from systemic. We ex- 
pected the features to have linewidths of ~ 20 30 km s -1 
FWHM, therefore a velocity resolution of 5-10 km s _1 
would have been sufficient for this project. To make the 
dataset generally usable to the community, we chose to 
observe at the full spectral resolution, with no on-line 
Hanning smoothing applied, and using the single IF cor- 
relator mode (sensitive only to right-handed circular po- 
larization). This resulted in 256 channels with a channel 
separation of 12.2 kHz (2.58 km s^ 1 ) and a channel width 
of 14.6 kHz (3.10 km s" 1 ). 

The uv calibrator data were edited in a standard it- 
erative way, with the inner 75% of the bandpass com- 
bined and inspected for obviously bad samples. Once 
all anomalous samples were excised from the calibrator 
data, the frequency-independent calibration was applied 
to the object data, which were then edited. The cali- 
bration solutions had a maximum closure error of 1% in 
amplitude and 1° in phase. The (frequency-dependent) 
bandpass correction was applied using the secondary cal- 
ibrator 3C 283 as reference, again following an iterative 
editing/calibrating method. Application of the bandpass 
correction produced closure errors of 0.5% in amplitude 
and 0.5° in phase, however because of the limited integra- 
tion time of the calibrator, the noise in the final data cube 
was increased by a factor of 1.8. With the editing and cal- 
ibration completed, the uv data were Hanning smoothed 
in velocity space. Besides improving the signal-to-noise 
in individual channels, smoothing with a Hanning ker- 
nel reduced the number of channels from 256 to 128 and 
made them independent, with channel width and sepa- 
ration equal at a value of 24.4 kHz (5.17 km s _1 ). 

The ultimate success of any interferometry observation 
hinges on whether a reliable map of source intensity in 
real space can be produced from incomplete sampling of 
the coherence function in the uv plane. At the imaging 
stage, weighting visibilities based on their local popu- 
lation density or distance from the tracking center af- 
fects the strength of sidelobe structure. We applied two 
weighting schemes to improve the quality of the prin- 
cipal s o lution : first, the "robust" scheme developed by 
Briggs] (|1995j ) to produce a mixture of natural and uni- 
form weighting with a single parameter; and second, 
a Gaussian taper function to downweight visibilities in 
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TABLE 3 

VLA Configuration and Data Parameters 



Dates of observation 
Time on source 
No. antennas 
Central freq. 
Primary beam HPBW 
Synth, beam HPBW 
Synth, beam PA 

1-<t noise (per 5.17 km s — 1 chan.) 
1-cr HI col. density (per 5.17 km s" 
1 mjy/beam 



chan.) 



24 Feb, 28 Feb, 1 Mar 1999 
12.2 hr 
26 

1.418 GHz 
31.'2 
37'.' 2 x 33'.' 1 
-72.5° 
0.94 mjy/beam 
4.6X10 18 cm~ 2 
0.49 K 



Raw data Processed/cleaned 



Bandwidth 

Velocity range 
No. channels 
Channel separation 
Channel width 



3.125 MHz 
660 km s 
183-843 km s" 1 
256 

2.58 km s- 1 
3.10 km s _1 



2.56 MHz 
540 km s _1 
230-770 km s" 1 
105 

5.17 km s" 1 
5.17 km s _1 



the sparsel y-sampled oute r uv plane. We used the MPS 
IMAGR task (jGreiscn 1998) and applied a variety of values 
for these two weighting schemes, empirically determining 
the best compromise between sidelobe structure, resolu- 
tion and sensitivity. The best compromise among all the 
considerations was a ROBUST value of and a Gaussian 
uw-taper with r. 3 = 6 kA, where r. 3 is the half-width at 
30% of maximum. This weighting scheme increased the 
noise in the final image by a factor of 1.2 compared to 
pure natural weighting and no icu-taper. 

A dirty image cube was produced from the full uv 
dataset, and channels free of H I line emission were iden- 
tified. The continuum was subtracted by fitting a first- 
order polynomial to 10 line- free channels on each edge 
of the band and subtracting this from all channels. A 
continuum-free dirty image cube was produced from this 
uv dataset, using a cell size of 7" per pixel and only in- 
cluding the central 105 channels to eliminate edge noise. 

To correct for incomplete sampling in the uv plane and 
reduce the level of the sidelobes, one typically uses the 
dirty image to iteratively construct a model of the true 
brightness distributi on. This is the strategy beh ind the 
CLEAN algorithm (jHogboml [l97l lOarkl Tl980h . under 
which the brightness distribution is modeled iteratively 
as a set of <5-function "clean components". CLEAN is 
very effective at removing sidelobes from point sources, 
since the clean components themselves are modeled as 
points. Spatially-extended H I emission, however, is 
poorly modeled by a set of point sources. Furthermore, 
while the strength of extended emission is usually rela- 
tively weak compared to a continuum point source, and 
thus the sidelobes associated with a single grid point of- 
ten fall below the noise in a dirty map, the summation 
of sidelobes due to complex extended structure can be 
significant and difficult to disentangle. 

To combat this problem, anoth er method of im- 
age de convolution was developed bv lWakker fc Schwarzl 
(1988). Their Multi-Resolution CLEAN (MRC) runs the 
CLEAN algorithm, but operates on the point-like and 
extended sources separately. First, the dirty map and 
beam are smoothed then subtracted from the original 
map and beam to obtain a full-resolution difference map 



and beam. Clean beams are constructed for the smooth 
and difference beams, and CLEAN is performed on the 
smooth and difference maps separately. By performing 
CLEAN on a smoothed map, two major improvements 
are made in reproducing extended emission: the clean 
beam is closer in size to the extended sources one is at- 
tempting to model; and the signal-to-noise is improved, 
allowing detection of fainter sources. The original res- 
olution is retained (and point sources are deconvolved) 
by cleaning the difference map separately and then com- 
bining the two resulting cleaned maps along with the 
full-resolution residuals. A detailed explanation of the 
method a nd the exact scaling param eters used are pre- 
sented bv lWakker fc Schwarzl (|1988fl . 

After trying standard CLEAN and MRC and compar- 
ing the results, it was determined that MRC was more 
effective in removing sidelobes. Our procedure was to 
clean the full region of each channel containing emission 
above the 3-c level, then examine the map for remaining 
strong sidelobes, which were present in most of the chan- 
nels. The output cube was spatially smoothed to twice 
the beam size and binned up by 5 channels, and regions 
that contained significant emission were manually delin- 
eated with polygons. The cleaning was repeated with 
MRC, using the original dirty cube and allowing clean 
components only within the specified region of each chan- 
nel. Successively deeper cleaning was performed, down 
to a threshold of 1-cr. 

The 1-tr noise level, calculated from the emission-free 
regions of each image plane, and other parameters for the 
cleaned data cube are summarized in Table [3] The noise 
of 0.94 mjy/beam is larger than the theoretical value 
of 0.7 mjy/beam for this instrumental setup and image 
weighting parameters, and this is likely a result of the 
bandpass calibration. The noise is similar from channel 
to channel, varying no more than 5%, as is shown in 
Figure [TJ 

To reduce image size and ease further computational 
analysis, the image cube was binned 2 x 2 in space. With 
a final pixel size of 14", the ~ 35" (FWHM) synthesized 
beam is still well-sampled. Noise was not reduced, how- 
ever, since neighboring pixels are correlated. The channel 
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Fig. 1. — The 1-a noise (crosses) and the peak H I flux (circles) 
in individual channels of the cleaned M 83 data cube. The noise 
levels have been multiplied by a factor of 100 and placed on the 
same scale as the peak flux to facilitate comparison. The noise is 
close to constant over the entire cube, with an RMS of 1.7% and 
maximum variations less than 5% of the mean. 

maps for the cleaned data cube are shown in Figure [5] 

2.2. Optical Data 

The optical data for M 83 were taken with the 0.6m 
Curtis Schmidt telescope located at Cerro Tololo Inter- 
American Observatory 5 (CTIO) in March 1999. Use of a 
thinned SITe 2K CCD chip resulted in a field of view of 
nearly 70' with 2.3" square pixels. To reduce flat-fielding 
errors that could be produced by non-repeatable filter 
wheel positioning and dislodged dust particles, we used 
a single Sloan r' filter, keeping the filter wheel unchanged 
throughout the run. A standard observing scheme was 
followed, with twenty bias frames and six to eight twi- 
light flats taken each evening, and standard star fields 
imaged at a range of airmass throughout the night. Ap- 
proximately half of the telescope time was spent imaging 
galaxies and half imaging blank sky flats, with exposure 
times of 15 minutes each. The total exposure times were 
10 hours on source and 8.75 hours on blank sky. Dark 
exposures taken during one cloudy night indicated no 
measurable CCD dark current. 

The initial reduction procedure followed a standard 
formula. The bias and overscan were removed from all 
twilight, sky and object frames. A first-order flatfield 
was created for each night from the twilights flats. This 
flatfield was applied to the sky flats, which were co- added 
and smoothed with a 128 x 128 boxcar to remove stars. 
This sky correction image was normalized and multiplied 
by the first-order twilight flat field to produce a master 
flat field for each night, which was applied to the M 83 
data. The individual galaxy frames were examined by 
eye, and eight of the 40 images contained undesirable 
features such as strong scattered light and were rejected. 
The remaining frames were sky-subtracted, aligned, and 
co-added. For proper noise estimation in future analysis, 
an average sky value of 4423 ADU was added back in to 
the co-added frame. Finally, the poorly-sampled outer 
regions were clipped, leaving a final 72'x69' image with 
a total integration time of 8.0 hours. 

Photometric calibration was accomplished with obser- 
vations of standard stars. The number of standard stars 
for the Sloan filter system was limited at the time of ob- 
servations, but using a star list provided by J. A. Smith 
(private communication) , we identified a number of stan- 

5 CTIO is part of the National Optical Astronomy Observatory, 
which is operated by the Association of Universities for Research in 
Astronomy (AURA), Inc., under cooperative agreement with the 
National Science Foundation. 



dard stars within lLandoltj (|1992f) fields. Nightly imaging 
of these fields yielded 50 observations of seven standard 
stars covering a range of airmasses (1.07 < secz < 2.0) 
and colors (0.0 < B — V < 1.0). For a pixel size of 
5.31 square arcsec and an effective exposure time of 900 
s, the surface brightness corresponding to 1 ADU pix -1 
in the final averaged galaxy frame was 30.11 ± 0.05 r' 
magnitudes per square arcsec. The average sky value of 
4423 ADU corresponds to a surface brightness of 21.04 
r' magnitudes per square arcsec, and the pixel-to-pixel 
noise in the final image (6.5 ADU) corresponds to 28.1 
r' magnitudes per square arcsec. Our 3-cr limiting sur- 
face brightness is 26.9 r' magnitudes per square arcsec, 
not accounting for systematic errors caused by scattered 
light and large-scale flatfielding effects. 

Foreground stars were removed by modeling the stellar 
point spread function (PSF) and m asking each star. Us- 
ing the DAOPHOT (jStetsonl [l987h package in IRAF 6 , 
we estimated the brightness of all stars in the frame 
by searching for 3-er peaks in the image and performing 
simple aperture photometry on them. This procedure 
produced 38,721 stellar objects; of these, 11 bright, iso- 
lated stars were chosen to construct the PSF. The PSF 
stars had an average (standardized) brightness of r'— 
16.3, with one moderately saturated star included to im- 
prove the fit in the wings. The PSF was constructed 
from the unsaturated regions of the PSF stars, using an 
elliptical Moffat function with (3 — 2.5 as the analytic 
PSF component plus a single look-up table for the resid- 
ual PSF component. The PSF was scaled to each PSF 
star and subtracted, leaving a residual image and reveal- 
ing previously-unidentified stars in the PSF stars' wings. 
These were fit with the PSF and removed from the orig- 
inal frame, and a new PSF was constructed. This proce- 
dure was iterated several times until the companion stars 
were completely removed and the PSF did not improve 
between iterations. 

Once the PSF was modeled, it was used to identify and 
measure all stars in the field. Of the stars originally iden- 
tified by aperture photometry, a subset of 33,950 were 
found to have profiles similar to the PSF. This amounts 
to an average of 6.8 field stars per square arcmin, down to 
a detection limit of r'= 23.8. The core of a star becomes 
saturated at r' ~ 14.5, but by fitting to the unsaturated 
wings, stars down to r' ~ 11 could be reliably measured. 

Stars were removed by masking each one out to the 
radius at which its flux fell below a specified threshold. 
This method is a compromise between elimination of pol- 
luting starlight and retention of usable data. To con- 
struct the mask, the model PSF was averaged radially 
to a 1-d function, scaled to each star according to its 
magnitude, and used to determine the radius at which 
the stellar flux fell below threshold. A threshold of 2 
ADU, well below the 7 ADU RMS noise in the co-added 
frame, was a satisfactory compromise between starlight 
removal and pixel retention. Typical PSF reference stars 
were masked out to a radius of 7 pixels (16"), and the 
brightest unsaturated stars were masked out to a radius 
of 9 pixels (21"), about 26 times larger than the average 

e IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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seeing HWHM of about 0.8". 

Some contaminating features did not conform to the 
shape of the PSF and were ignored by the above proce- 
dure. These included background galaxies, very bright 
saturated stars, diffraction spikes, and charge transfer 
("bleeding") features, and they were masked by hand. 

3. THE NORMAL GALACTIC DISK OF M 83 
3.1. H I Moment Map Construction 

The bulk distribution and kinematics of emission-line 
gas are most easily summarized by moment maps, in- 
cluding the total H I surface brightness (0th moment) 
and the intensity-weighted mean velocity field (1st mo- 
ment). To construct the moment maps, the image cube 
was first smoothed spatially to a round beam of 52" 
FWHM, roughly 1.5 times the size of the original synthe- 
sized beam. This smoothing is done to improve the S /N 
and ensure that the faint extended components of emis- 
sion peaks are included in the summation. Pixels with 
flux greater than 4-cr were included, as were immediate 
neighbor pixels in both space and velocity. 

The moment maps are shown in Figures [3] and 0J The 
primary beam attenuation correction was applied to only 
the moment map (which is shown in column den- 
sity units) after integrating in velocity. Throughout this 
work, all H I intensity maps displaying either integrated 
H I emission or H I column density have been corrected 
for the primary beam sensitivity. 

The total H I mass determined from the H I sur- 
face brightness map is 4.7xl0 9 Mq. This is lower 
than the distance-corrected va l ue of 6 .1 x 10 9 obtained by 
iHuchtmeier fc Bohnenstengell (fl981 ). who mapped the 
emission out to about twice the radius as the present 
study. Their single-dish scans were impervious to zero- 
spacing issues of interferometry and recovered all of the 
emission in the beam. Our result is therefore in reason- 
able agreement with their value. 

3.2. HI Disk Features 

The total H I map (Figure [3]) is a good tracer of spa- 
tial morphology, and it shows a number of notable fea- 
tures. The inner disk of M 83 contains a minimum in 
the H I approximately 1' (1.3 kpc) in diameter, with a 
measured column density of 4xl0 20 cm~ 2 . This feature 
is surrounded by the inner gaseous disk, which extends 
to a radius of 6' (8 kpc) and ranges in column density 
from a minimum of 7xl0 20 cm~ 2 to a (probably opti- 
cally thick) maximum of 1.7xl0 21 cm -2 . The inner disk 
shows a hint of spiral structure, with an irregular pattern 
of bright spots, and an aspect ratio close to unity, as one 
would expect for a nearly face-on galaxy. The small-scale 
structure is unresolved in our map, but it agrees quali- 
tatively with the results of iTilanus fc Allenl |l993), who 
studied M 83 at higher spatial resolution. 

Beyond 5', the H I brightness falls off abruptly, with 
A/hi ;$ 4xl0 20 cm~ 2 in the region immediately around 
the inner disk. In the regions NE and SW of the inner 
disk, the column density becomes very low, approaching 
5xl0 19 cm~ 2 . Further out, there is evidence of a bright 
ring of H I at a radius of 10' (13 kpc). The ring has a 
peak column density of lxlO 21 cm~ 2 , similar to the in- 
ner disk, and the projected thickness varies between 60" 
and 80". With the smoothed beam size (52") this scale 



is unresolved, but the width can be no more than 1' (1.3 
kpc) after removal of the instrumental profile. The emis- 
sion in the ring appears clumpy, with the brightest region 
in an arc on the west side of the galaxy. In addition, the 
aspect ratio of the ring is quite different from that of the 
inner disk, with semi-major and semi-minor axes of 10' 
and 5' respectively. The position angle is close to that 
of the inner disk (see below) . 

Outside of the ring, the H I spatial distribution be- 
comes irregular, with arms and streams at various pro- 
jected orientations. One arm appears just beyond the 
ring on the NW side, and spans about 25' (33 kpc) to 
a point 16' (21 kpc) NE of the disk center. There is a 
corresponding arm on the SE to S side of the galaxy, at a 
similar projected distance. Two notable features are the 
arms further west of this, which extend from a common 
point 14' (18 kpc) west of the center and extend toward 
the north. The longer of these arms reaches a projected 
galactic radius of 26' (34 kpc) and is clearly detached 
from the rest of the disk. Both NW arms appear to be 
mirrored in the SE, and similar mirroring can be seen 
between other features and clumps in the extended H I. 
Since the primary antenna sensitivity falls to 10% of the 
central value for the most extended arms, these features 
could be even more extended than observed. It is clear 
that the arms are distinct from the main H I disk, down 
to our detection limit. The symmetrical nature and vary- 
ing orientations of the outer arms indicates either a warp 
that varies with radius or a collection of kincmatically 
distinct orbiting streams. 

Analysis of the kinematic morphology can improve our 
understanding of the complex outer disk features, and 
this is most thoroughly done by combining the velocity 
field map, the individual channel maps, and the data 
cube as a whole. The velocity map is presented in both 
Figure[3] (as contours) and Figured (as an image map and 
contours). In both figures, the inner disk shows a typi- 
cal differential rotation pattern (a so-called "spider dia- 
gram"), with a central systemic value of ~ 500 km s . 
To estimate the systemic velocity and orientation, we fit 
the central 5' with a Brand t rotation curve using the GAL 
task in MPS (|Greisenlll998l) . We find a systemic velocity 
of «hci = 513 ± 2 km s _1 , where the error is estimated 
from the spread of best-fit values at different radii, and 
using the optical definition i>hci,opt = c (A bs — A re f)/A re f. 
The best-fit position angle is 225° ± 2° (measured to the 
receding side from north to east), and the best-fit incli- 
nation is i = 23° ± 2°. 

Our estimates of the inner H I disk parameters are 
consistent with previous results. Previous values for the 
system i c velo city range from 505 km s" 1 determined by 
IComtei (|1981[) through Ha interf e romet ry to 516 km s _1 
quoted by Ide Vaucouleurs et al.1 (|199lD and determined 
by H I observations. The HIPASS Brig ht Galaxy Cata- 
log r eports a value of 513 ± 2 km s^ 1 (Koribalsk i et al.l 
2004). We adopt 513 km s" 1 as the systemic velocity for 
the remainder of this work. Our estimates of the orienta- 
tio n and inclination are co nsistent with the those quoted 
by ITilanus fc Alien] (|1993[ ). who have data with higher 
spatial resolution. Therefore we adopt i = 22° for our 
inclination, and note that the maximum rotation rate of 
the inner disk is u ro t,max = (whci.max - w sys )/ sini = 200 
km s -1 . 
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Fig. 2. — M 83 H I channel maps. Heliocentric velocities are marked in the upper right of each channel map. The grayscale ranges from 
—3 to 15 mjy to show the fainter structure in each channel, therefore most of the bright features appear saturated. Only channels with 
measurable H I emission are displayed. 
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Figure [2] (continued) 
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Outside of the inner disk, a discontinuity develops in 
the rotational velocity. This is most clearly seen in Fig- 
ure [3l where the velocity contours bend abruptly into a 
warp, and in the channel maps for 440 km s _1 < fhciio < 
540 km s^ 1 (see Figure [5]). On the W side of the galaxy, 
the warp correlates well with the H I ring in projection, 
but on the E side the warp peaks inside of the ring. The 
isovelocity contours are shifted by 2' counter-clockwise 
from where they would be under ordinary differential 
galactic rotation. Beyond the ring, the isovelocity con- 
tours switch direction again, and they no longer reflect 
a flat rotation curve. In addition, the aspect ratio of the 
disk increases in this region, indicating that the axis of 
rotation is becoming more inclined to our line of sight 
as the radius increases. The highest radial velocities are 
seen near the very edge of the detected H I distribution, 
at a distance of 8' (11 kpc) N and S of the galactic center. 

It is apparent from this analysis that the inner disk, 
the ring, and the outer arms are kinematically distinct 
features. The ring and outer arms could be consid- 
ered anomalous features themselves, however we wish 
to search for anomalous- velocity emission superposed on 
the bulk H I emission, and for this purpose we assume 
the outer arms are part of an extended disk of varying 
inclination. 

3.3. Large-Scale Optical Features 

The final optical image for M 83 shows many of the 
typical features of a face-on, grand-design spiral galaxy 
(see Figured]). The spiral arms are seen radiating from 
a central bar-like structure, and these arms contain a 
number of star clusters and dust lanes. The integrated 
H I intensity traces the spiral arms closely. The spiral 
structure is apparent only out to a radius of about 4.5' 
(6 kpc) in the north and 5.5' (7 kpc) in the south, at 
which point the surface brightness has dropped to 22.5 



r 1 magnitudes per square arcsec. This is the same radius 
at which the inner H I disk falls off abruptly and at which 
the warp in the velocity field appears. 

Beyond the inner optical disk, there is an extended 
optical envelope devoid of any structure. This envelope 
reaches a surface brightness of 26 r' mag per square arc- 
sec at a radius of about 10' (13 kpc), falling to unde- 
tectable levels at the location of the H I ring. In addi- 
tion, the aspect ratio of the extended light, appearing 
stretched northeast to southwest, is different from that 
of the inner disk, which is nearly circular. 

An arc-like feature can be seen projected 20' (26 
kpc) northwest of center, extending 12' (16 kpc) north- 
east to southwest. This feature is barely visible 
in the Palomar Sky Survey plate covering this re- 
gion, and it has been cataloged as a dwarf galaxy by 
iKarachentseva fc Karachentsevl (|1998f ). who designate it 
KK208. The surface brightness across the feature varies 
from 25-26 r' mag per square arcsec, although there are 
three bright stars along the arc that complicate masking 
and flux measurements. There is no measurable H I any- 
where along the feature down to a 3-c column density 
limit of 1.5xl0 19 cm~ 2 . The arc lies between two outer 
H I arms, and it has a shape similar to those arms. 

4. ANOMALOUS- VELOCITY H I: EXTENDED 
EMISSION 

4.1. HI Disk Modeling 

To search for anomalous-velocity material, one must 
first establish the criteria that make such material 
"anomalous" . The difficulty of this can be appreci- 
ated by considering the ensemble of Galactic HVCs. 
A commonly-used definition limits membership to com- 
plexes with velocities deviating by greater than 100 
km s _1 (or sometimes 90 km s _1 ) relative to the Local 
Standard of Rest (LSR). This definition was intended to 



10 



Miller et al. 




separate the HVCs from the so-called Intermediate Ve- 
locity Clouds (IVCs) and gas in the plane. However, such 
a definition does not take into account differential galac- 
tic rotation, an effect which can vary the observed veloc- 
ity of a cloud by as much as 100 km s" 1 , depending on 
the line of sight direction (|Wakkerj|1991l ). To account for 
this effect, some authors have used the velocity relative 
to the so-called Galactic Standard of Rest (GSR), which 
is a point at the location of the Sun that does not par- 
ticipate in galactic rotation. This definition works well 
for gas at high \z\. However, it over-corrects the velocity 
of any cloud that has a rotational compo nent of motion , 
such as may be found closer to the plane. Wakkcr (1991) 
introduced the concept of "deviation velocity," which is 
the difference between the observed LSR velocity of the 
cloud and the maximum velocity expected along that line 
of sight due to differential rotation, assuming a simple 
model for the gaseous Galactic disk. At low latitudes in 
the Milky Way this varies greatly, but at latitudes above 
~ 40° the velocity limit is around 60 km s _1 . 

For external spiral galaxies, and especially face-on 
galaxies, the situation is simplified in that the bulk disk 
velocity is easily measured along any line of sight. Mate- 
rial at velocities different from this, assuming some veloc- 
ity threshold or kinematic model for the disk, would be 



kincmatically distinct from the cold H I disk and there- 
fore anomalous. Hereafter we define anomalous-velocity 
material to be any H I emission that does not partici- 
pate in the normal differential rotation of the bulk of H I 
along that line of sight. Implicit in the definition is the 
inclusion of spatially isolated H I clumps, since the disk 
rotation is undefined in regions where the disk does not 
exist. Therefore an anomalous- velocity clump (hereafter 
AVC or AV gas) can be anomalous in velocity or position. 

As low- velocity AV gas could overlap and merge with 
the disk velocity, it is useful to remove any H I disk 
emission. Additionally, it is important to ignore this 
emission in a non-interactive source detection scheme. 
We modeled the disk of M 83 by fitting a Gaussian 
to the spectrum along every pixel of the cleaned im- 
age cube. The fitting was performed only along lines 
of sight where the peak exceeded 3.5-cr (3.3 mJy/beam), 
and only included channels whose flux exceeded 1.7 -a 
(1.6 mJy/beam). These criteria were applied to ensure 
that faint, velocity-broadened emission would not con- 
taminate the thin disk fit. 

Disk fitting and removal failed at warp discontinuities, 
where the H I disk has multiple velocity components 
across a scale of about 1'. It is not clear whether this 
spatial overlap of kinematically distinct disk components 
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is a real feature or an effect of beam smearing, but in ei- 
ther case it poses a difficulty for effective disk removal. 
It is clear that the warp features are not AV material 
in our definition, since they define the disk velocity. We 
therefore attempted to remove this emission by refitting 
the original data cube with multiple Gaussian compo- 
nents. For lines-of-sight where the residual flux matched 
the above single-Gaussian criteria (i.e., peak above 3.5- 
cr, 3 or more adjacent pixels with flux above 1.7-a), a 
two-component Gaussian was fit. The fitting threshold 
was kept high so that true AV gas would escape the fit. 

To estimate the quality of the disk model, we exam- 
ined the reduced \ 2 an( i covariance matrices of all the 
fits performed. The fits that included only a few data 
points were obviously poorly constrained, but inspection 
of the residual data cube showed that even these fits were 
reasonable. There were several lines-of-sight that had no- 
ticeably bad fits, but these were generally isolated pixels 
whose neighbors had similar parameters to each other. 
To remove the outliers, we performed a median box fil- 
tering on the fitting parameter images, using a kernel 
width of 3 pixels. The disk model and residual cube 
were reproduced from the adjusted fitting parameters, 
and inspection showed that the poor fits were removed 
and the disk subtraction was improved. 

The fitted line-of-sight velocity dispersion varies across 



TABLE 4 

Extended AV Gas: Dynamical Properties 



Component 


F 




KE 




(Jy km s~ 1 ) 


(10 7 M Q ) 


(10 53 erg) 


Full residual cube 


203.0 


97.0 


21.0 


Inner disk 


17.6 


8.4 


16.0 


Inner disk, low-u 


11.8 


5.6 


8.7 



the galaxy, as can be seen in Figure [H It reaches a 
maximum of 20 km s _1 in the center, averages about 
15 km s _1 in the central H I disk and ring, and falls to 
10 km s _1 or less in the extended disk and arms. 

4.2. Extended, Disk-Like AV Gas 

Removal of the H I disk results in a residual cube with 
significant H I emission remaining. For ease of display, 
we created position-velocity (jp-v) plots, or slices through 
velocity along a line of specified location, length, and po- 
sition angle. The slices were taken parallel to the major 
axis at a spacing of 70" (twice the beamwidth). The 
value of a given pixel in the p-v image equals the distance- 
weighted average of the nearest four pixels in the channel. 
A sample of slices is shown in Figure [7] for both the full 
data cube and the residual. 
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cm 2 and increasing by factors of 2. Also labeled is KK208, a faint companion galaxy or stellar stream to the 



Comparison of the slices with and without the disk 
emission clearly shows an emission component not well 
modeled by differential galactic rotation. This is es- 
pecially evident along the major axis. The anomalous 
component is spatially extended and clumpy, with devi- 
ation velocities generally ±40-50 km s _1 near the emis- 
sion peak along any line of sight. The emission gener- 
ally appears brighter on one velocity wing of the disk 
H I, toward the systemic galactic velocity 7 . Some of the 
anomalous emission occurs in the "forbidden" velocity 
region, which is the p-v quadrant on the opposite side of 
the systemic velocity from the disk. One interpretation 
is that material at velocities between the disk and galac- 
tic systemic would be rotating more slowly than the bulk 
of material in the galaxy. Material in the forbidden p-v 
region would be moving in the opposite sense of the bulk 
rotation along that line of sight, i.e., counter-rotating. 
These interpretations assume the disk and anomalous 
material have similar spatial distribution and rotation 
axis orientation. In a face-on galaxy like M 83, material 

7 Hereafter, velocities on the systemic side of the disk velocity 
will be referred to as "low relative velocity", while velocities on 
the opposite side of the disk velocity will be referred to as "high 
relative velocity". This designation places the velocities in the 
reference frame of the rotating galaxy. 



in the forbidden region is likely to be in vertical motion. 

To quantify the amount of AV gas in M 83, we con- 
structed column density and mean velocity maps from 
the residual cube (see Figure [8]) . The dynamical prop- 
erties of the AV H I were determined from these maps 
and are summarized in Table 2] The mass of AV H I in 
the residual cube is 9.7x10 s M Q , corresponding to 20% 
of the total H I. It was apparent that poor subtraction 
of the H I ring/warp and outer H I arms contaminated 
the residual cube, therefore we manually defined a re- 
gion containing only anomalous emission projected on 
the central disk. The region used was an ellipse centered 
on the galaxy with a semi-major axis of 6.0' (radius of 
8 kpc), a semi-minor axis of 4.5', and a position angle 
of 45° as shown in Figures [5] This area contains the 
most dynamically homogeneous disk gas and produced 
reasonable model parameters. The total H I mass in this 
region, including the cold disk, is 1.02xl0 9 Mq, of which 
8.4xl0 7 Mq or 8.2% is present in the anomalous compo- 
nent ("inner disk" in Tabled]). The low-relative- velocity 
emission contributes 5.6xl0 7 Mq or 5.5% of the total 
inner disk H I. 

The velocity field shows two high- velocity spots to the 
east and west of the galaxy center, identified by rect- 
angles in Figured] These features are also seen on the 
high relative velocity side of the p-v slices, including one 
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Fig. 6. — Line-of-sight H I velocity dispersion for M 83, obtained by fitting a Gaussian model to the disk emission. The fit was performed 
along every line of sight (pixel), but the values shown have been median box filtered to remove poor fits. The velocity dispersion is 20 
km s _1 in the center, falling to less than 10 km s _1 in the outer regions. 
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-20 -10 10 20 -20 -10 10 20 
arcmin arcmin 

Fig. 7. — Position- velocity (p-v) diagrams for M 83, taken parallel to the major axis at separations of 70", as plotted on the H I emission 
map (upper right) with the slice direction indicated by the arrow. Each slice, boxcar smoothed in velocity by 3 channels, is shown with the 
disk emission included (top) and subtracted (bottom). The rotation curve is traced by the white and black points. Grayscale ranges from 
to 6 mjy/beam. Black and white contours begin at +2-cr and increase by factors of two. Red contours begin at — 2-cr and decrease by 
factors of two. Residual emission can be seen in all the slices, especially between the systemic and rotational velocity along the major axis. 
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structure +1' along the +140" p-v diagram in Figure [3 
at velocities between +50 and +100 km s _1 . In a simple 
model of disk rotation, these clumps would be rotating 
more quickly than the bulk of the H I. They are discrete 
in projection compared to the low relative velocity com- 
ponent described above, and will be discussed in detail 
in Section E] (as AVCs 1 and 2). 

Aside from these features, the motion and morphology 
of the residual emission mimic that of the galactic disk, 
with the southwest side receding and the northeast side 
approaching. The line-of-sight velocity dispersion ranges 
from 10-15 km s _1 , similar to the thin disk, although this 
may be an underestimate if the wings of the anomalous 
emission are excluded in the thin disk mask. A map 
of deviation velocity was produced by subtracting the 
velocity map of the modeled disk from that of the residual 
cube (see Figure [5]), and it shows a rotation rate 40-50 
km s _1 slower in projection. If the inclination of this 
disk of AV emission is also about 20°, then its rotation 
rate is about 100-150 km s -1 , compared to 200 km s _1 
for the cold disk. 

We used the deviation velocities to construct a map of 
the kinetic energy with respect to the cold disk (see Fig- 
ure [5]). Unlike the previous analysis, this interpretation 
assumes all motion is perpendicular to the galaxy plane 8 , 
and it provides an upper limit to the energy apportioned 
to vertical bulk motion. The total kinetic energy so de- 
rived is 2.1 xlO 54 erg, with 1.6 xlO 54 erg from gas within 
the inner disk and 8.7xl0 53 erg of that from the low- 
relative-velocity material. The last figure excludes dis- 
crete clumps, such as the one projected 2' northwest of 
center; these features possess nearly half of the kinetic 
energy of the inner AV gas, and are fully discussed in Sec- 
tion[5] The 1.6xl0 54 erg from the inner 8 kpc is equiva- 
lent to the energy of 8 supernovae kpc -2 , although the 
z component of velocity is not known for this material, 
so this is an upper limit to the kinetic energy. 

While the distribution of extended gas is well- 
characterized by moment maps and p-v slices, faint, com- 
pact emission sources are easily overlooked. Detection of 
such systems requires a quantitative, objective searching 
technique that utilizes the entire set of data in a statisti- 
cal way. The next section addresses our search for faint, 
discrete H I emission sources. 

5. ANOMALOUS- VELOCITY H I: DISCRETE 
EMISSION 

5.1. Searching for AVCs 

Compared to other observational subfields of astron- 
omy, radio synthesis imaging suffers from a lack of robust, 
objective source detection software. Identifying sources 
and measuring their fluxes and other parameters is typi- 
cally done by hand from moment maps and similar data 
products. For bright, point-like sources, this is not a 
problem, but difficulties arise when one wishes to search 
statistically for emission of unknown scale and location 
that may be projected on much stronger sources. 

We have developed a suite of software that searches a 

8 Lacking any knowledge of transverse motion, we actually as- 
sumed motion fully along the line of sight, which is offset by 24° 
from the galactic pole direction. Thus velocities and kinetic en- 
ergies are 10% and 20% lower, respectively, then what would be 
calculated for gas with only ^-direction motion. 



3-d data cube for signal in a statistical way. The soft- 
ware is described in detail in Appendix [A] and summa- 
rized here. The data cube is smoothed with a tunable 
velocity filter to emphasize sources of a particular ve - 
locity width, following the method of lUson et alJ (|1991h . 
Peaks in this smoothed cube are identified and used to 
define 3-d "islands" of emission by way of contouring, 
using a variation of the CLUMPFIND algorithm devel- 
oped bv lWilliams et ail (|1994f ). The velocity smoothing 
allows detection of kinematically broad but statistically 
significant features that might otherwise be overlooked. 
In addition, this algorithm allows a quantitative analysis 
of errors with the use of simulated datasets. 

To search for discrete sources, we masked the residual 
H I cube within ±20 km s -1 of the H I disk velocity to 
remove poorly-subtracted emission in the center of the 
disk profile. We ran the software with a peak thresh- 
old of 4-er, a clump merging threshold of 2.5-cr, and a 
clump extending threshold of 1-a. The width of the 
Gaussian filtering kernel was varied from 3-9 channels 
(15-45 km s _1 ). For sources that were detected at mul- 
tiple filter widths, we used the results from the filtering 
that produced the highest S/N (i.e., the kernel that was 
closest in width to the unsmoothed feature). The spa- 
tial and kinematic morphology of the candidate clouds 
was analyzed interactively using a variety of tools, in- 
cluding p-v plots, channel maps and 3-d visual rendering 
tools. A two-component Gaussian fit along the velocity 
axis at the clump centroid recovered the velocity and line 
width of the clump, as well as the line width of any H I 
disk emission along that line of sight. Finally, the flux of 
each clump was verified by hand using the standard radio 
method of aperture photometry. Values were consistent 
with the results from our automated software. 

5.2. Detected Systems 

We discovered 14 discrete clouds of H I emission at 
anomalous velocities or positions in the M 83 data cube. 
These anomalous-velocity clumps (AVCs) are distinct 
from the bulk of the H I disk emission in that their 
contours of emission detach from the disk at levels ap- 
proaching the map noise. The AVC properties are given 
in Tables [5] and [H The locations of the AVCs are shown 
superposed on the H I column density map in Figure [10] 
and superposed on the velocity map in Figure 1111 

The AVCs are found across the H I disk of M 83, but 
they do exhibit a small degree of clustering, and several 
are located in regions where the kinematics of the disk 
are most complicated. Three clumps (1, 2 and 4) are lo- 
cated within the kinematically "normal" inner disk. One 
clump (3) is found at the inner disk/warp interface, two 
additional sources (7 and 8) are projected just beyond 
the ring, and four clumps (5, 6, 9, and 12) are found pro- 
jected near outer H I arms. The remaining four AVCs 
(10, 11, 13, 14) are located well away from the center of 
the galaxy, at projected radii where the primary beam 
sensitivity falls to 10% or lower. 

Gaussian fits and p-v plots were produced as described 
in the previous section, and these are shown for each 
AVC in Figure [T3J The integrated H I intensity maps 
for the spatially-extended AVCs are shown in Figure [T3l 
To eliminate the confusion of spatially- varying noise, the 
p-v slices and intensity maps have not been corrected for 
the primary beam attenuation. The fluxes, masses, and 
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Fig. 8. — Maps of the AV emission in M 83, showing (left) column density (in units of 10 20 cm -2 ) and (right) mean velocity. The 
ellipse indicates the inner, well-modeled portion of the cold H I disk, and the rectangles identify discreet high-relative- velocity features, as 
described in the text. The anomalous gas appears disk-like, rotating about 40-50 km s _1 more slowly in projection than the cold H I disk. 




aremm arcmin 

Fig. 9. — Maps of the AV emission in M 83, showing (left) deviation velocity and (right) kinetic energy per 14" X 14" pixel (in units of 
10 50 erg). The contours show the mean H I disk velocity in steps of 10 km s — 1 for reference, and the other notations are as described in 
Figure[8] The anomalous gas appears disk-like, rotating about 40—50 km s — 1 more slowly in projection than the cold H I disk, but with 
about the same rotation axis. Most of the kinetic energy is contained in the feature 2' northwest of center. 
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TABLE 5 

Discrete AVCs: General Properties 



AVC a 


RA 


Dec 


w AVC b 


n c 






Vr 


.e 




(J2000) 


(J2000) 


(km/s) 


(sq. arcmin) 


(') 


(kpc) 


(') 


(kpc) 


1 


13 36 49.2 


-29 51 13 


595± 2 


2.63 


5.3 


6.9 


1.7 


2.2 


2 


13 37 19.3 


-29 51 13 


398± 2 


2.75 


5.5 


7.2 


3.7 


4.8 


3 


13 36 42.7 


-29 57 31 


667± 2 


0.77 


2.8 


3.6 


7.1 


9.3 


4 


13 37 18.3 


-29 49 35 


547± 2 


0.58 


2.5 


3.3 


4.4 


5.7 


5 


13 37 12.9 


-30 04 59 


512± 2 


0.39 


1.7 


2.2 


13.4 


17.4 


6 


13 37 08.6 


-30 06 23 


498± 3 


0.39 


1.8 


2.3 


14.6 


19.0 


7 


13 37 06.4 


-30 02 11 


411± 2 


0.39 


1.5 


2.0 


10.0 


13.0 


8 


13 36 56.8 


-29 43 31 


361± 2 


0.39 


1.9 


2.5 


8.4 


11.0 


9 


13 37 33.4 


-30 10 35 


690± 8 


0.39 


1.3 


1.7 


19.9 


25.8 


10 


13 38 11.8 


-29 27 38 


593± 5 


0.39 


1.2 


1.6 


28.6 


37.2 


11 


13 37 21.6 


-30 15 15 


407±12 


0.39 


1.6 


2.1 


23.9 


31.1 


12 


13 35 45.7 


-29 50 58 


591± 5 


0.39 


1.0 


1.3 


16.5 


21.5 


13 


13 37 25.9 


-30 17 35 


405± 5 


0.39 


1.1 


1.8 


25.9 


33.7 


11 


13 35 19.0 


-29 37 39 


633±32 


0.39 


1.2 


1.6 


26.3 


34.2 



a AVCs 1—8 arc considered real detections, AVCs 9-14 arc likely spurious. 

Heliocentric radial velocity of the peak of the AVC emission. 
c Solid angle subtended by the AVC after collapsing its HI emission along the velocity axis. A 
value of 0.39 sq. arcmin indicates the feature is unresolved. 

Greatest linear extent of the AVC, uncorrected for broadening by the synthesized beamwidth. 
c Projected distance from the galactic center. 



TABLE 6 

Discrete AVCs: Dynamical Properties 



AVC a 


TP b 
-fincas 


r? c 

-tcorr 


M m d 




fdcv f 




KE h 




(Jy km/s) 


(Jy km/s) 


(10 6 M Q ) 


(10 20 cm" 2 ) 


(km/s) 


(km/s) 


(10 53 erg) 


1 


3.08 ±0.07 


3.12 ±0.08 


14.9 ±0.4 


2.9 


+51 


29±2 


3.9 ±0.1 


2 


1.42 ±0.06 


1.48 ±0.06 


7.1 ±0.3 


1.1 


-41 


28±2 


1.2 ±0.0 


3 


0.99 ±0.05 


1.12 ±0.05 


5.4 ±0.3 


2.6 


+74 


56±4 


2.9 ±0.1 


4 


0.23 ±0.03 


0.24 ±0.03 


1.2 ±0.1 


0.6 


+113 


32±2 


1.5 ±0.2 


5 


0.11 ±0.02 


0.18 ±0.04 


0.8 ±0.2 


1.0 


-77 


15±2 


0.5 ±0.1 


6 


0.11 ±0.02 


0.21 ±0.04 


1.0 ±0.2 


1.1 


-103 


15±2 


1.0 ±0.2 


7 


0.13 ±0.02 


0.17 ±0.03 


0.8 ±0.2 


1.2 


-166 


44±3 


2.3 ±0.4 


8 


0.11 ±0.02 


0.13 ±0.03 


0.6 ±0.1 


0.7 


-90 


24±2 


0.5 ±0.1 


9 


0.08 ±0.02 


0.24 ±0.06 


1.2 ±0.3 


2.3 




30±2 




10 


0.14 ±0.03 


2.77 ±0.52 


13.2 ±2.5 


29.1 




24±2 




11 


0.19 ±0.03 


1.09 ±0.17 


5.2 ±0.8 


9.2 




82±5 




12 


0.06 ±0.02 


0.14 ±0.04 


0.7 ±0.2 


1.3 


+135 


29±2 


1.2 ±0.3 


13 


0.16 ±0.03 


1.51 ±0.27 


7.2 ±1.3 


13.5 




20±2 




14 


0.15 ±0.03 


1.68 ±0.31 


8.0 ±1.5 


17.4 




164±12 





a AVCs 1—8 arc considered real detections, AVCs 9-14 arc likely spurious. 

Total flux in the AVC, uncorrected for the primary beam attenuation. Errors arc 1-cr. 
c Total flux in the AVC, corrected for the primary beam attenuation. 
d Mass of HI contained in the AVC. 

G Peak column density of HI, integrated across the velocity width of the emission. 

^ Deviation velocity of the AVC, defined as the difference between the radial velocity of the peak AVC emission 
and the velocity fitted to the HI disk (if detectable) at that location, as described in Section |5.1| 
s Full width at half maximum of a Gaussian fit to the velocity profile. 

Kinetic energy of the AVC found from the deviation velocity and mass. This is a lower limit to the total kinetic 
energy as only radial velocities are used. 
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Fig. 10. — The H I column density map for M 83, showing the locations of the detected AVCs. The shapes indicate the projected spatial 
extent of each AVC out to l-cr in H I column density (5xl0 18 cm -2 ). The large gray circles indicate the primary beam sensitivity. The 
column density is in units of 10 20 cm -2 . 

column densities reported in Table [5] have been corrected 
for this effect. 

A number of the AVCs have fluxes or velocity structure 
in the p-v plots that suggest they are spurious. From a 
"by-eye" analysis, we estimate that eight of these features 
(AVCs 1-8) are real detections and the remainder are 
not. Below we describe specific properties of individual 
AVCs, grouping them according to similar distribution 
and attributes. For completeness, we include features 
thought to be spurious, although these are excluded in 
our discussion in Section [6l 
AVC 1 

This source is the brightest and most massive AVC de- 
tected, with M = 1.5xl0 7 Mq. Centered 2 kpc from the 
dynamical center of M 83, AVC 1 is also the only discrete 
clump projected over the inner H I disk. The intensity 
map and p-v plots show that it is spatially resolved and 
extended in velocity, with a deviation velocity Udcv = 51 
km s _1 and line width of 29 km s" 1 FWHM. The mea- 
sured line width is likely only accurate near the emission 
peak, as the p-v slices display a connection in velocity 
to the underlying disk emission. The projected intensity 
map supports this idea, with a secondary maximum pro- 
jecting to the south, where the cold disk approaches the 



velocity of the cloud. The characteristics of this AVC 
are consistent with material flowing into or out of the 
disk, although the lack of spatial information along the 
line of sight complicates conclusions about the source's 
z height. 

The moment and dynamical maps described in Section 
14.21 clearly show emission from this feature (see Figures 
[HlandEi]). The peak of the emission can be seen 1.5' west 
and 1.5' north of center in Figure [8^,. 
AVCs 2, 4 

These two clumps are close together in projection at a ra- 
dius of about 5 kpc (near the edge of the inner disk) , yet 
they differ greatly in their other characteristics. AVC 
2 is bright, massive (M = 7xl0 6 Mq) and extended 
(fi = 2.8 sq. arcmin, I — 7 kpc), while AVC 4 is less 
massive (M = lxlO 6 Mq) and barely resolved (largest 
linear extent I — 3.3 kpc). In addition, the deviation ve- 
locities of these clumps are quite different, with AVC 2 
at high relative velocity (vdev = —41 km s -1 ) and AVC 
4 at low relative velocity (t>dcv = 113 km s" 1 ). The line 
widths are similar, but the emission of AVC 2 appears 
to blend with the disk velocity. Emission from AVC 4 
also blends kinematically with the disk, but the peak is 
more clearly separated, as can be seen in the velocity 
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Fig. 11. — The H I intensity- weighted mean velocity map for M 83, showing the locations of the detected AVCs. Each square is centered 
at the centroid of the AVC, with the color indicating its heliocentric radial velocity. 



profiles in Figure [12] The heliocentric velocity of AVC 
4 (522 km s _1 ) places it in the "forbidden" region of 
counter-rotating material in the p-v plane. The spatial 
coincidence of these clumps hints at a correlation be- 
twee n them, possibly in the form of an expanding bubble 
(e.g.. lKamphuis et al.lll99lh . but it is difficult to reconcile 
their differing morphological and kinematic structure. As 
with AVC 1, these two clumps appear in the AV gas maps 
shown in Figures [8] and [9] 
AVC 3 

At a projected distance of 9 kpc from the galactic cen- 
ter, this AVC appears in a region where the H I ring 
is disturbed and has a relatively low column density of 
A/hi = 3xf0 20 cm~ 2 . The clump is extended in space, 
covering two synthesized beams in area and spanning 3.6 
kpc. It is also extended in velocity, with a line width of 
56 km s _1 FWHM, and its emission appears to merge 
with the disk H I at lower velocities. The source is mas- 
sive (M = 5.4 xlO 6 Mq) and has a large kinetic energy 
of 3xl0 53 erg. 

Although the clump is projected on a region of chang- 
ing inclination between the inner disk and outer arms, it 
does not appear to result from beam-smearing of a warp 
discontinuity. First, it is a compact structure, whereas 
a warp feature would be extended like the H I ring or 



outer arms. Second, it is much broader in velocity than 
the kinematically cold H I elsewhere in the galaxy. Fi- 
nally, the kinematics of the surrounding disk material are 
regular in comparison to that seen at warp interfaces in 
other regions of the galaxy. This cloud is the only such 
feature seen along the H I ring. 

This AVC appears at a location of low disk Ahi- The 
p-v plots and Figure [10] show a lack of emission in the 
disk a few arcmin to the east of this location, although 
H I disk emission is present at the projected location of 
AVC 3. 
AVCs 5, 6 

These two AVCs have similar mass, size, line width, and 
deviation velocity. At a projected radius of about 14' 
(18 kpc) to the south of the galactic center, AVCs 5 and 
6 are separated by 2' (2.6 kpc) in projected space and 
20 km s _1 in velocity. AVC 6 is the more massive and 
kinematically anomalous clump, with a mass of lxlO 6 
Mq and deviation velocity of —103 km s _1 . AVC 5 has a 
mass of 8xl0 5 M Q and deviation velocity of —77 km s _1 . 
The line width of both clumps is 15 km s _1 , producing 
a velocity dispersion of 6 km s _1 . Unlike the previously 
discussed AVCs, the emission from these H I clumps does 
not merge with that of the disk in position or velocity. 
In addition to similar kinematics, these two clouds have 
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Fig. 12. — Position-velocity slices taken along the E-W direction and velocity profiles are shown for each AVC, which is centered at the 
intersection of the solid lines in the p-v plot. Grayscale ranges from to 6 mjy/beam. Black and white contours begin at 2-<r (cr = 0.94 
mjy/beam) and increase by factors of two. Red contours begin at — 2-cr and decrease by factors of two. The profile plot, taken at the 
centroid of the AVC total H I emission, shows the Gaussian velocity fit to the AVC and H I disk (if present). AVCs which show up in 
multiple slices are identified by small numerals. 
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Figure 1121 (continued) 
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Fig. 13. — Integrated H I intensity maps for the four spatially- 
extended AVCs. AVCs 1 and 2 have bright cores and low-level 
extended features, while AVCs 3 and 4 are more compact and only 
mildly resolved. Grayscale ranges from 0-0.25 Jy beam -1 km s -1 , 
and the contours are spaced at 0.025 Jy beam -1 km s — 1 . The 
ellipse in the AVC 3 plot shows the HPBW of the synthesized 
beam. Coordinates are with respect to the center of M 83, with 
north up and east to the left. 

similar sizes and are spatially unresolved by the VLA 
synthesized beam. The FWHM of a Gaussian fit across 
either feature is 45"; when deconvolved from the ~ 35" 
synthesized beam, this implies a maximum source size of 
30" (~ 0.7 kpc). Their similar characteristics and loca- 
tion hint at a possible relation, and they might be bright 
spots in a single H I feature, although this speculation is 
made less likely by the 20 km s _1 difference in deviation 
velocity and the small radial velocity dispersions. 
AVCs 7, 8 

These emission sources are similar in that they appear 
between the H I ring and an outer arm, although they 
are projected on opposite sides of center. Both have 
highly negative deviation velocities (vdev = —166 and 
—90 km s _1 ), are spatially unresolved, and have similar 
fluxes and masses (M = 8 and 6xl0 5 M Q ). Clump 7 has 
a large fraction of the kinetic energy of the AVC ensem- 
ble, due primarily to its large deviation velocity. From 
the integrated H I maps, which have the same colormap 
scale, one can see that AVC 7 is more sharply peaked 
than AVC 8. Consideration of the p-v plots and velocity 
profiles shows that the line of sight to AVC 7 contains 
emission at a variety of other velocities (e.g., —100, —40, 
+40 km s _1 ). The clumps producing this emission fall 
below our threshold for clump identification, and it is 
unclear whether they result from real velocity structure 
in the AVC or from a systematic problem such as poor 
continuum subtraction or bandpass calibration. 
AVCs 9, 12 

Two clumps lie superposed on or near H I arms outside 
the half power width of the primary beam. Both are 
spatially unresolved and have similar line widths of 30 
km s" 1 FWHM. At these radii, the primary beam at- 
tenuation correction factor is about 3, so although both 
have small measured fluxes, the corrected masses are sig- 
nificant (1.1 and 0.6xl0 6 M G ). Both AVCs have high 



TABLE 7 
Optical Counterparts to HI AVCs 



AVC a 


b 

/V,o 


/V,llm C 


Notes 


1 


20.4 




inner disk, spiral arm 


2 


20.7 




inner disk, spiral arm 


3 


24.2 




inner disk 


4 


22.2 




inner disk 


5 


27.8 


27.2 


stars 


6 


27.9 


27.4 


stars 


7 


25.4 


27.1 


stars 


8 


24.4 


26.4 


outer disk, stars 


9 


26.5 


27.0 


outer disk, stars 


10 


26.8 


27.3 


stars, galaxy 1.5' SE 


11 


27.2 


26.8 


stars 


12 


27.5 


26.8 


bright coincident stars 


13 


28.1 


27.4 


stars 


14 


27.7 


26.9 


stars, galaxy 



a AVCs 1—8 arc considered real detections, AVCs 9-14 are 
likely spurious. 

The median surface brightness, in r' magnitudes per 
square arcscc, at the location of the AVC. For extended 
AVCs, the median was taken within the HPBW centered 
on the peak HI emission. 

c The 3-cr surface brightness limit on any optical features 
coincident with the AVC. This was determined by fitting 
and subtracting a plane from the optical image at that 
point. See that text for details. 

The notes identify optical features coincident with the 
AVC, including foreground stars, extended emission from 
M S3 (designated by location in the disk), and background 
galaxies. 

deviation velocities, with Vdev = +135 km s _1 for clump 
12, and a velocity difference of +100 km s _1 for clump 
9 measured from the H I arm projected nearby. These 
detections are the least statistically significant of the set 
(4-ct and 3-cr for AVC 9 and 12, respectively), and they 
are hereafter treated as spurious. 
AVCs 10, 11, 13, 14 

All four of these clumps are in low-sensitivity regions, 
where the attenuation correction factor approaches 10. 
AVCs 11 and 13 are projected close together in space 
(2' apart) and velocity (2 km s^ 1 apart). All four have 
strongly peaked emission maps, and their lines of sight 
contain emission at velocities throughout the data cube. 
In addition, AVCs 11 and 14 have very large line widths 
compared to the other AVCs. It is likely that these AVCs 
are not real, but are the result of systematic effects such 
as calibration errors. If they are real, the beam correc- 
tion produces large masses and kinetic energies for these 
clumps. 

5.3. Optical Counterparts to the AVCs 

Searching for optical counterparts to the H I AVCs was 
carried out by hand. Each AVC was overlaid on the un- 
masked optical image, which was inspected for optical 
emission (see Figure fl4|) . To search for small-scale emis- 
sion structure superposed on the extended galaxy profile, 
the masked optical image was median filtered with a 5 x 5 
pixel box. The location of each AVC was fit with a plane, 
excluding the four AVCs within the inner disk, which is 
difficult to fit with a simple model. The fit residuals 
were inspected for signal above the noise, which was not 
found in any case. The noise in the residuals was used 
to estimate surface brightness upper limits for optical 
counterparts. The results are shown in Table [71 



24 



Miller et al. 



O 
crj 



20 



10 







-10 



-20 



-30 



N. 



OlO- 



- ~E- 




3 XP 




O 3 



Cm' r*"" t/' , * Lr 



*W i J 1 i - 



■ .. - - s ;^-^ wl-Q kp c'w 
j I L_j i i I i i [ i i i lJ l12 i i I i i l_j I 



- - 26 



28 



24 



20 



10 







10 



20 



30 



arcmin 

Fig. 14. — The optical M 83 image overlaid with the locations of the H I AVCs. 
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Below we describe possible counterparts and limits to 
the optical surface brightness for the individual AVCs. 
We discuss these limits in the context of other deep op- 
tical studies in Section lrT2"1 
AVCs 1, 2, 3, 4, 7, 8 

These AVCS are all projected over the bright optical disk. 
The bright cores of AVCs 1 and 2 are projected on spiral 
arms, appearing to trace the shape of the arms, although 
this is possibly due to the method used to detect and de- 
lineate the H I clumps. The optical light coincident with 
AVC 1 peaks at 20.4 r' mag per square arcsec, exclud- 
ing the bright star clusters in the spiral arms, and falls 
off to 22.1 r' mag per square arcsec in the outer part of 
the clump. Likewise, the core of AVC 2 is projected on 
a bright spiral arm, with ji i — 20.7 r' mag per square 
arcsec. The optical light falls to 24.5 r' mag per square 
arcsec in the outer regions of this extended H I clump. 
While these AVCs align with spiral arms that contain 
star clusters, the spatial resolution of the H I data is in- 
sufficient to draw direct comparisons between AVC and 
star cluster morphology. 

The other two AVCs within the inner disk, AVCs 3 
and 4, have no obvious small-scale optical counterparts 
aside from stars. The inner edge of AVC 4 coincides with 
the outer edge of a spiral arm, although this is again 
likely due to the searching and masking technique we 
have employed. AVC 3 is projected on a disruption in 



the H I disk, but no evidence of such a feature is seen in 
the optical data. 

Clumps 7 and 8 are projected further out, over the 
south and north edge of the extended stellar envelope. 
The underlying surface brightnesses in these regions is 
25.4 and 24.4 r' mag per square arcsec, respectively. 
There are no small-scale optical features, aside from 
stars, down to the detection limit of about 27 r' mag 
per square arcsec. 
AVCs 5, 6, 9, 11, 13 

The majority of the large-radius AVCs contain only stars 
within the lowest H I contour. Clumps 5 and 6 are iso- 
lated from bright stars, but AVCs 9, 11 and 13 lie in the 
wings of very bright stars. As a result, much of the light 
in these regions is masked out to determine the surface 
brightness, and the amount of flux below a limiting mag- 
nitude of about 26 r mag per square arcsec is poorly 
constrained. There is no evidence for emission down to 
this level, however. 
AVCs 10, 12, 14 

These three AVCs are positioned on or near optical fea- 
tures that could be related. AVC 10 is projected 1.5' 
northwest of an object with a galactic light profile. The 
contours of AVC 14 contain a similar object. Both ob- 
jects are smaller than the extent of the AVCs, although 
the H I clumps are unresolved. Given their apparent 
size, it is likely that these optical sources are background 
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galaxies, although it is impossible to verify without spec- 
troscopic data. 

Clump 12 is centered on the location of two bright 
stars, making the limiting magnitude at this location dif- 
ficult to determine. There is no signal at this position in 
the 1.4 GHz continuum map, therefore it is unlikely that 
the detected H I results from poor continuum subtrac- 
tion. It is more likely that this is a chance superposition 
between bright stars and a detected H I source (real or 
spurious) . 

5.4. AVC Detection Simulations 

To analyze our results quantitatively, it is necessary 
to determine significance of the detections in terms of 
the completeness and false detection rate. Interferometry 
maps differ from images made with non-interferometric 
devices in that neighboring spatial pixels are not inde- 
pendent. The transformation of limited samples from 
the uv plane to the image plane results in the signal 
from every location being spatially smeared. Resolution 
elements therefore overlap each other, and more impor- 
tantly, noise features are extended and often similar in 
size to the objects being searched for. Prior to primary 
beam correction, the noise is nearly constant across the 
field. It is upon this background of spatially-dependent 
fluctuations that one searches for discrete emission fea- 
tures. 

To investigate the sensitivity of our data, we simulated 
2000 observations of Gaussian noise in the uv plane, us- 
ing the same uv distribution and noise characteristics 
as our filtered M 83 observations. Each dataset was 
transformed to the image plane using the same technique 
described in Section 12.11 so that the resultant image 
cube had dimensions, spatial coherence (i.e., synthesized 
beamwidth) and a identical to that of the cleaned, binned 
M 83 data cube. Each simulated cube was processed with 
our detection software, using identical parameters as de- 
scribed in Section IBTTl except for a single filtering kernel 
of 5 channels (25 km s" 1 ) FWHM, the average linewidth 
of the AVCs. The histogram of clump flux values, scaled 
to the parameter space volume of a single masked M 83 
data cube, is shown in Figure [T5j The low-flux side of 
the histogram shows a rapid fall-off due to decreasing 
detection efficiency. The histogram peaks near 0.045 Jy 
km s _1 , which is about 4.5 times the l-er flux for a single 
beam in the smoothed cube (see Appendix [A} ■ This is 
consistent with our detection threshold of 4-c 

We can use the simulated distribution to calculate the 
likelihood of obtaining at least one spurious detection 
above a given flux threshold in the data cube. If we 
assume the false detection rate is determined by Pois- 
son statistics, we find that the likelihood of obtaining 
at least one detection above flux F is given by pf = 
1 - e~ tJlF n F /v\ = 1 - er^ F for v = 0. The mean ex- 
pected number of detections /i is given by the integral 
of the empirical distribution between F and +oo, and it 
is a straightforward matter to find F for various values 
of pf- The flux limits for pp = 0.32, 0.05, and 0.01 are 
plotted in Figure \T5\ Several of our AVCs have fluxes 
(uncorrected for attenuation) within the region of non- 
zero false detection rate. If we exclude the three signif- 
icant clumps with uncorrected fluxes near or above 1.0 
Jy km s _1 (AVCs 1-3), we are left with 11 clumps in the 
0.06-0.23 Jy km s _1 flux range. These are indicated in 
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Fig. 15. — Flux distribution of false detections in the simulated 
pure noise data cubes. The binned values have been scaled to equal 
the mean number of false detections per masked M 83 data cube 
per 0.01 Jy km s flux bin, so that the integral of this distribution 
equals the mean expected number of false detections in our data 
cube, a value of 4.1. The dotted lines note the likelihood pp of 
obtaining at least one false detection above the given flux value in 
a single data cube (see text for further explanation). Open circles 
indicate the fluxes of detected clumps in the M 83 data, with the 
three brightest clumps (AVCs 1-3) excluded. For clump flux > 0.1 
Jy km s - , there is an excess over purely spurious detections. 
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Fig. 16. — Detection efficiency of the simulated clumps. The 
efficiency is 90% above 0.15 Jy km s _1 . 

Figure [T5l 

Integrating the full empirical distribution leads to an 
expected mean of 4.1 false detections per masked data 
cube. With an algorithm flux detection threshold of 4-cr, 
and for 177,500 independent samples in our masked data 
cube, we expect 5.3 detections of greater than 4-cr signif- 
icance from purely statistical considerations. That the 
simulations average fewer false detections indicate limi- 
tations of our detection algorithm. Both false detection 
rates are consistent with our "by-eye" analysis indicating 
that AVCs 9-14 are spurious. The chance of obtaining at 
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least 6 spurious detections given the simulation mean of 
4.1 is 22%; given the statistical mean of 5.3, the chance 
is 44%. 

We note that the simulations do not account for sys- 
tematic errors which may occur from bandpass calibra- 
tion or continuum subtraction and which may introduce 
artifacts similar to clumps of emission. Some of our "de- 
tections" show indications of this effect (e.g., the multi- 
ple velocity components of AVCs 11 and 14). Excluding 
these, the number of false detections is still consistent 
with the expected rate. Modeling of systematic calibra- 
tion errors will be included in future versions of the sim- 
ulation software. 

The completeness of our sample was determined from 
the detection efficiency, which in turn was found from ad- 
ditional simulations. Twenty data cubes with noise char- 
acteristics similar to the data were created in the manner 
described above. To each of them were added 200 3-d 
Gaussian sources with spatial FWHM = 40" (approxi- 
mately the synthesized beamwidth), velocity FWHM = 
5 channels (25 km s _1 ), and peak brightness varying be- 
tween 1.0 and 5.0 mJy/beam in steps of 0.5 mJy/beam. 
The integrated fluxes of the sources varied between 0.01 
and 0.27 Jy km s _1 . The positions of the sources were 
produced semi-randomly so that the centers were not 
necessarily at integer pixels but the sources did not over- 
lap within 5-er in any direction. The resulting cubes 
were processed with our detection software in an iden- 
tical fashion as the source data. 

At least 90% of the simulated sources were detected to 
a flux of 0.16 Jy km s _1 , as is shown in Figure \W[ Be- 
low this, the detection efficiency falls off rapidly, reaching 
50% at a flux of about 0.1 Jy km s _1 . The shape of the 
efficiency curve may explain the lack of low-flux detec- 
tions in our sample. Also from these simulations, we note 
a systematic flux overestimate of ~ 10% for a simulated 
0.2 Jy km s _1 clump and ~ 30% for 0.1 Jy km s _1 . The 
measured AVC fluxes have not been corrected for this 
effect. 

6. DISCUSSION 

Different HVC production scenarios predict different 
characteristics for the HVCs. The galactic fountain pre- 
dicts that material will be projected on the disk, while 
tidal stripping and accretion predict HVCs across the 
field of view. We can divide the detected AVCs into two 
categories: those projected on the disk (AVCs 1, 2 and 
4 and the low relative velocity extended emission), and 
those projected away from the disk. Separate treatment 
of each group places constraints on the importance of the 
various HVC formation schemes. 

6.1. Galactic Fountain: "The Beard" and HVCs 

Under the galactic fountain model, hot gas is ejected 
from the disk by multiple supernova explosions in star- 
forming regions. This material rises at roughly the adi- 
abatic sound speed, moves outward due to a decrease 
in the radial gravitational potential, and rotates more 
slowly due to conservation of angular momentum. After 
the gas cools, it loses its buoyancy an d spirals back to 
the d i sk, either as sm all clouds or sheets (jShapiro fe Fieldl 
[19761 iBre gmanlll980l ). 

The observational implications of this model depend 
in detail on the nature of the hot corona and mag- 



netic field, which are ill-constrained, but several gen- 
eral features should be apparent in observations of ex- 
ternal galaxies. First of all, discrete clouds of scale ~ 
1 kpc should be observable at velocities between sys- 
temic and systemic plus the coronal sound speed (< 150 
km s _1 ), on both wings of the disk velocity profile. Sec- 
ond, H I emission could appear at velocities between 
systemic and the disk, i.e., on the low-relative velocity 
wing of the disk profile. This emission could be spatially 
extended (spanning a few kpc) if small-scale condensa- 
tions are suppressed (lField|ll965t iMathews fc Bregmanl 
[19781: iChevalier fc Oegerlei 119791 iBregmanl (T980f T FT 
nally, anomalous velocity H I should coincide with the 
spiral arms and Ha disk emission, especially the discrete 
high- velocity clumps, which would be infalling near their 
point of origin. 

Our results are consistent with these predictions. The 
extended anomalous disk of H I appears to be rotating 
more slowly than the main H I disk, which would be the 
case for emission from a galactic fountain. This mate- 
rial must be spatially separated from the bulk H I, and 
we conclude that it lies in a vertically extended disk. A 
similar phenomenon has been observed in a number of in- 
clined spiral galaxies, inclu ding NGC 89 1 (ISwaters et al 
19971 lOoster loo et al.l l200"l. NGC 2403 (IFraternali et al. 
2002b]), NGC 4559 jBarbieri et all l2005h. NGC 253 
(iBoomsma et al.ll2005ft . and NGC 6946 (jBoomsma et al.l 
|2008j) . In particular for the the i = 60° spiral NGC 
2403, IFraternali etafl (|2002bh argue that this emission, 
termed the "beard" due to its appearance in p-v dia- 
grams, arises from a vertically extended H I component 
that is rotating 20-50 km s _1 more slowly than the cold 
disk. They model the anomalous disk and characterize 
it as having a radial inflow at velocity < 20 km s _1 , im- 
plying an inflowing mass rate of M ~ 0.3-0.6 Mq yr" 1 
for an infall timescale of ~ 10 8 yr. They also detect 10 6 - 
10 7 Mq "streams" and "spurs" of H I emission, ranging 
in size from 5-10 kpc and deviation velocity from 50-60 
km s _1 . They conclude that this material arises from 
a galactic fountain, with a mass exchange rate close to 
the predicted value of ~ 1 Mq yr -1 . Further evidence 
is provided by Chandra X-ray observations showing an 
extended corona of T ~ 5xl0 6 K gas with a radiativ e 
cooling rate of 0.1-0.2 Mr* vr" 1 (Fraternali eT al. 2002a). 
Additional modeling bv IFraternali fc Binnevl (|2006l ) sup- 
ports this conclusion, although it suggests that some 
amount of IGM accretion must be present to temper the 
fountain and produce the low angular momentum mate- 
rial. 

The face-on inclination of M 83 makes it difficult to 
model t he rotation of the beard c o mponent in the same 
way a s IFraternali et al.l (|2002bh . IFraternali fc Binnevl 
(2006), and others. It is also not clear whether the 
anomalous disk is radially infalling as is the case of NGC 
2403, nor is it clear how vertically extended the material 
is. If the beard is the result of a galactic fountain, it 
is likely that it represents material near maximum z be- 
cause of its low velocity, and it probably is gas that has 
recently condensed and is beginning to fall back onto the 
disk. The discrete clouds projected on the disk at higher 
deviation velocities would be in freefall at lower z. If the 
column density of a cloud exceeds that of the disk, the 
cloud would punch through and disrupt the H I disk in 
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that location. AVC 3 is projected over a region of low 
disk Nm, however it is projected near the warp and not 
the spiral arms, so it may have a different origin than the 
other disk-coincident AVCs. 

The dynamical properties of the AV gas are consistent 
with what would be expected in a galactic fountain. The 
beard mass of 5.6xl0 7 M Q (8.4x10 7 M q including the 
AVCs) and an estimated cooling/freefall time of 5xl0 7 
yr leads to a mass exchange rate of 1 Mq yr" 1 , similar 
to the predictions of the model (|Bregmanl [T9801 . The 
star formation rates of M 83 and the Milky Way are sim- 
ilar, so one w ould expect similar mass exchange rates. 
iTalbotl (|1980h estimates a star formation rate (SFR) in 
the Milky Way between 1.7 M@ yr -1 (from converting 
CO measurement s to H 2) and 2.0 Mq yr -1 (from Ha). 
iBell fc Kennicuttl (|2001f ) constrain the the SFR of M 83 
to be between 1.1-2.4 Mq yr" 1 from F UV and Ha data, 
using the relations of iKennicuttl (|1998h . The mass infall 
rates in both galaxies are estimated to be ~ 1 M Q yr -1 , 
consistent with the similar SFR. The kinetic energy of 
the discrete AVCs projected over the disk total about 
7xl0 53 erg, equaling the total output of about 1000 su- 
pernovae. This is a lower limit to the kinetic energy 
initially imparted to the gas by the collection of super- 
novae, since the velocity is measured in only one direction 
and the z height of the AVCs is unknown. A supernova 
rate of 0.01 yr -1 is consistent with this if only 1% of the 
energy is converted to kinetic energy. 

6.2. Galactic Disruption and Accretion: External 
HVCs 

It is difficult to reconcile the presence of AVCs pro- 
jected outside the stellar disk with the galactic fountain 
model, thus there is likely another HVC production sce- 
nario at work. AVCs 3, 7 and 8 are projected on the 
H I warp/ring, and are possibly associated with that. 
The remainder of the AVCs are compact, and all except 
AVCs 5 and 6 are likely spurious detections. 

The optical results rule out emission from a population 
of dwarf galaxies similar to those in the Local Group. 
The 3-cr detection limits for the AVCs projected outside 
the H I disk are at least 26.8 r' magnitudes per square 
arcsec. Assuming an integrate d color V — R ~ 0.4-0.6, 
typical of Local Group dwarfs (Mateo 1998), an d usin g 
the SDSS filter transformation from lSmith et all (2002), 
this corresponds to about 27.2 V magnitudes per square 
arcsec. This limit is faint enou gh to detect th e center of 
any of the dwarfs tabulated bv lMateol jl998). Likewise, 
the inferred ratios of Mm/iy > 2M Q /L© are larger than 
those of the majority of Local Group dwarfs, indicating 
a dearth of starlight compared to these systems. 

O ur optical depth is sim ilar to that of other HVC stud- 
ies. iWillman et all <|2002) use SDSS data to search for 
stellar components in 13 Galactic HVCs, detecting none 
to limits of 26.7-3 0.1 V magnitudes per square arcsec. 
ISimon et all (|2006f) obtain an upper limit of 25.25-26.25 
V magnitudes per square arcsec for HVC Complex H, 
depending on the assumed stellar population. Fainter 
companions and stellar streams are observed in the Milky 
Way and M 3 1, but they do no t have evidence of H I 
emission (e.g- IZucker et al.ll2004t IWillman et ail 12005) . 

If they are real, these H I clouds are best explained 
by a combination of tidal disruption and accretion. The 
H I warp, H I ring, and elongated, gas-free companion 



KK208 are all evidence supporting a tidal interaction at 
some point in the galaxy's past. In addition, an unpub- 
lished wide-field H I map shows evidence of an H I arm 
to the north of the galaxy, extending east fr om the tip 
of th e longest NW arm visible in Figure [3] (jPark et al.l 
2001). Our field of view is insufficient to include this 
region, but it is further evidence of a tidal interaction 
that can strip diffuse material and produce anomalous 
H I emission. 

6.3. Limits on the Mass Distribution of Galactic HVCs 

Our results are sensitive enough to place constraints on 
the mass distribution of the HVC ensemble in M 83 and, 
by extension, the Milky Way. Because of the variable 
mass sensitivity across the field, we restrict our analy- 
sis to AVCs detected within the HPBW of the primary 
beam. We assume that the full distribution of M 83 
HVCs falls within this region, which is a reasonable as- 
sumption if the vertical distance z is less than about 
50 kpc (and zsini < 17 kpc or 15'). Only the eight 
confirmed detections lie within the primary HPBW, and 
we scale our estimates of the false detection rate to the 
reduced number of independent sightlines we are now 
sampling. Five of the eight detected clumps have masses 
in a region where the false detection rate is significant 
(greater than 0.1 detection expected per mass bin), while 
the remaining AVCs are much more massive and un- 
likely to be spurious. Assuming a Poisson distribution 
for false detections, and using the integrated false de- 
tection rate as the mean "background" expected for a 
mass bin, we determine the background-corrected mass 
distribution with 90% confidence limits. This is shown 
in Figure [T7l using the tabulated data of lGehrelsl (1986) 
to calculate confidence limits for low number counts. We 
obtain upper limits above an HVC mass of 1.2xl0 6 M© 
and below a mass of 6xlO 5 M ; below 4xlO 5 M , the 
completeness falls below 50%, so we have cut off the dis- 
tribution at this point. 

To draw conclusions about the Galactic HVC popula- 
tion, two simplifying assumptions must be made. First, 
we assume that the underlying HVC mass distributions 
are identical for the Milky Way and M 83. Second, we 
assume the Galactic HVCs are contained in a uniform 
spherical distribution about the Galactic center. While 
this is likely incorrect, it is a sensible choice given our lack 
of knowledge of their distances. The collection of Galac- 
tic HVCs across the whole sk y has been cataloged by 
Wak ker fc van Woerdenl (|1991h . and excluding the Mag- 
ellanic Stream and outer Galactic arm, we randomly as- 
sign distances drawn from the assumed distribution to 
the cataloged HVCs. The assigned distances and cat- 
aloged H I fluxes produce an estimate of the Galactic 
HVC mass distribution. This procedure was repeated 
500 times to fill in the high- mass regions, and the result- 
ing distribution was scaled to the number of HVCs in the 
catalog (560). 

The derived HVC mass distribution is plotted in Figure 
[17] for four different mean Galactic HVC distances, along 
with the HVC mass distribution for M 83. While the 
confidence limits are large, it is clear that an ensemble 
distance of less than about 25 kpc is consistent with the 
M 83 distribution, while a larger average distance results 
in a larger number of massive clouds. Such clouds would 
have been easily visible in our M 83 observations. It is 
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Fig. 17. — Mass distribution of the Milky Way HVCs, assuming 
uniform spherical distributi ons of varying mean dista nce and using 
the HVC catalog of IWakker fc van Wocrdcn (1991). The points 
show the distribution of our lowest-mass detections in M 83, with 
95% confidence errorbars and upper limits. Three detections lie 
off the horizontal scale at Mpji > 5xl0 6 Mq. This distribution 
has been corrected for the false discovery rate (the "background" ) 
and the source detection efficiency. If the underlying HVC mass 
distributions of the two galaxies are identical, the Milky Way HVCs 
must be closer than about 25 kpc or we would have detected more 
HVCs in M 83. 



possible that some HVCs in M 83 could fall outside the 
HPBW and therefore be missed if the distance is much 
greater than 50 kpc, however even at 50 kpc we detect an 
order of magnitude fewer HVCs than should be present 
in M 83. In addition, the heliocentric distances are up- 
per limits to the more fundamental z height of the HVCs, 
therefore the distances of the HVCs from the plane are 
consistent with being le ss than about 10 kpc. Finally, i f 
the HVCs cataloged by IWakker fc van Woerdenl (| 19911 ) 
are not distinct but rather clumps within larger com- 
plexes, then the mass distribution derived here will un- 
derestimate the number of high-mass HVCs, and the con- 
straints imposed by our results would be stronger. 

These constr aints are consistent with previous results. 
IWakkerl (|2001h has published a catalog of 18 HVCs 
and 8 IVCs with distances and metallicities constrained 
by absorption against background stars or extragalactic 
sources. The metallicities vary from solar to a few per- 
cent of solar. The upper bounds to the vertical z heights 
range from 0.1 kpc to 7 kpc, and the corresponding mass 
limits rang e from about 10 4 to 2xl0 6 M©. More recent 
results indicate distance s of 5-15 kpc to a handful of ad- 
ditional Ga lactic HVCs (iThom et alj|2q06t IWakker et all 
120071 I2008D . Similarly. iPutman et all (|2003[) have de- 
tected Ha emission in 25 HVCs and use this to limit the 
distance of these objects from the ionizing Galactic ra- 
diation field and thus constrain their height above the 
Galactic plane. The values range from 5 < z < 40 kpc. 
Limits from external galaxies are similar. iThilker et al.l 
( 2004) detected an ensemble of 20 discrete H I clouds 
within 50 kpc of M 31. These sources have H I masses 
in the rang e 10 5 -10 7 M Q . A handful of the objects iden- 



tified by these authors would be visible if present in M 
83, and our results are consistent with this. iPisano et al.l 



(|2007f ) survey six analogs of the Local Group and detect 
no extragalactic H I emission. Using a similar analysis 
to ours, they conclude the Galactic HVCs must be dis- 
tributed within 90 kpc and have average H I mass less 
than about 4x 10 5 M Q . The handful of discrete HVCs we 
detect in M 83 could represent the high-mass tail of a 
much more abundant low-mass, nearby population. 

7. SUMMARY AND CONCLUSIONS 

We have performed deep, wide-field imaging of the H I 
in M 83, mapping the outer fea tures for the firs t time . 
The H I ring and warp seen by iTilanus fc Alle"n] (|1993l ) 
are confirmed, and the outer arms extend to a radius 
of 34 kpc and show reflective symmetry, hinting at a 
tidal formation scenario. At least 80% of the H I mass 
is located outside the optical extent of the galaxy, with 
the possibility that more H I lies undetected outside the 
rimary beam HPBW. The optical c ompanion KK208 
Karachentseva fc Karachent scv 1998, see Figure [SJ con- 
tains no H I to a limiting column density of 1 x 10 19 cm~ 2 . 

We discovered a spatially-extended component of 
anomalous- velocity gas deviating 40-50 km s _1 from the 
bulk H I disk and coincident with it in projection, with 
a line-of-sight velocity dispersion of 10-15 km s _1 . We 
interpreted this as a vertically extended disk rotating in 
the same sense but about 100 km s _1 more slowly than 
the kinematically cold, thin disk. The 5.6xl0 7 Mq of 
H I it contains is 5.5% of the total H I within the stellar 
disk of the galaxy. 

We have introduced a new technique of searching for 
faint extended radio emission, combining several exist- 
ing tools and including statistically robust data mod- 
eling. Other recent studies have used similar tech- 
niques of contourin g and extended sou rce extraction 
(jPutman et al.ll2002t lde~Heii et al.ll2002ft . These devel- 
opments are driven by the need for lower detection lim- 
its in projects constraining the missing satellites problem 
and obtaining a full census of the Galactic HVCs. 

Using this technique, our observations are sensitive to 
unresolved sources (d < 1 kpc) with masses greater than 
5xl0 5 Mq. We have discovered 14 discrete anomalous- 
velocity emission sources, of which we consider 8 to be 
real detections, ranging in H I mass from 6xl0 5 to 
1.5 x 10 7 Mq and projected on and off the disk. Three of 
these sources are high-significance detections and appear 
to be distinct H I clouds with masses in excess of 5xl0 6 
Mq. They are spatially extended and coincide in projec- 
tion with the optical spiral arms. The flux distribution 
of the low-mass sources is inconsistent with purely spuri- 
ous detections, therefore we conclude that most of these 
are real HVCs in the M 83 system. They are generally 
unresolved, off-disk features and are free of diffuse opti- 
cal light to a limiting surface brightness of 27 r' mag per 
square arcsec. 

We conclude that a combination of a galactic fountain 
and tidal stripping are responsible for the anomalous H I 
features that we observe in M 83. The handful of HVCs 
we have detected are likely to represent the high-mass 
end of a substantial HVC population in this galaxy. We 
expect future deep H I observations of nearby spirals to 
add to the growing number of HVC analogs detected in 
external galaxies, and to thereby shed light on the ubiq- 
uity of HVC activity and the nature of our own Galaxy's 
gaseous neighborhood. 



High- Velocity Clouds in M 83 



29 



We thank the staff of the VLA and CTIO for their ous referee for constructive comments that improved the 

assistance with the planning and execution of these manuscript. EDM would like to thank Mario Mateo, 

observations. We especially thank Jonathan Williams Hugh Aller, and Tim McKay for their helpful sugges- 

and Juan Uson for providing access to and help with tions toward the improvement of this work, 
their source detection software code, and the anonmy- 

APPENDIX 

A NEW SOURCE DETECTION METHOD FOR SPECTRAL SYNTHESIS IMAGING 

We have developed a suite of software that searches a spatially-cohere nt, 3-d data cube f or sig nal in a statistical way . 
This software borrows heavily from the methods and code developed bv lUson et all (jl991l ) and I Williams et alj (j!994t) . 
but as we have made variations to the methods and translated much of the code to a new programming language, we 
here describe the new software in some detail. 

The program begins by smoothing the input data cube in velocity with Gaussian kernels of the form 

H(v) = e -(" 2 /2^) ) (Al) 

where a v is the width of the particular Gaussian. The convolution is performed by a fast Fourier transform. If one 
assumes a Gaussian velocity profile for a clump of emission, the signal due to this feature as a function of channel or 
frequency v would be 

S(v) = A(v )e- (v ~ vo)2/2a *.° , (A2) 

and the smoothing function acts as a matched filter, enhancing features of similar width and de-emphasizing narrower 
or wider features. In a given channel, the normalized amplitude of the smoothed feature is given by 



S(v)H(v)dv 

(AM) = J — , (A3) 

H{vfdv 



where the integral is taken over the entire feature. This value is maximized when a v = <j v .o- The normalization factor 

1/ J H(v) 2 dv is such that the amplitude of the smoothed feature will equal the amplitude of the original feature under 

this same condition; otherwise, the smoothed amplitude will be depressed compared to the unsmoothed value. 
One can calculate the significance level of such a detection by recognizing that the smoothed noise is 



<^M 2 ) = 



<Ti(v) 2 H(v — vo) 2 dv 

72 . (A4) 

H{v) 2 dv 



if we assume that neighboring channels are independent and the noise (given as Oi for channel i) sums in quadrature. 
The noise in each channel is determined by reflecting the negative portion of the data histogram across the mean and 
performing an iterative 4-cr rejection until the measured a converges. The significance level (S/N) is calculated as 
(A(u ))/(cr(wo) 2 ) 1 ^ 2 and written to disk as a data cube. 

Peaks in S/N are identified down to some detection threshold, typically around 4-cr, with neighboring above-threshold 
pixels merged into groups. These 3-d groups or "islands" serve as kernels for an iterative contouring scheme, with 
successively lower thresholds applied to the S/N cube. Each group is extended by adding new pixels that are above 
threshold and share at least two corners with any member pixel. That is, only neighbors that share at least one 
dimensional plane with the member are added, so that neighbors at a diagonal in all th ree dimensions a re exc luded. 
The algorithm for this procedure is based on the CLUMPFIND software developed by I Williams et al.l (JT994) . The 
contours are produced in intervals (which only affect the processing speed) down to a specified level, typically 3-cr, 
iterating at each level until no new neighbors are added to any groups and merging groups which meet the neighbor 
criterion. Finally, any immediate neighbors (i.e., sharing four corners or two parameter planes) above 1.5-ct are added 
to each group to improve flux measurements. No merging or iterating is done at this point. Groups are eliminated if 
they contain fewer pixels than a completely unresolved source (10 pixels in our case), as these features are unphysical. 

The group IDs are written to a data cube, which is applied to the original data cube as a mask for each group in 
turn. The total flux, centroid, and other parameters of the group are determined from the original data cube, using 
only the pixels that are members of the group. The primary beam correction is applied at this point, and the corrected 
flux and mass are recorded along with the uncorrected parameters. 
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